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Abstract 

Recent data have determined that the structure of the high pressure e phase of solid oxygen 
consists of clusters composed of four O2 molecules. This finding has opened the question about 
the nature of the intermolecular interactions within the molecular oxygen tetramer. We use multi- 
configurational ab initio calculations to obtain an adequate characterization of the ground singlet 
state of (02)4 which is compatible with the non magnetic character of the e phase. In contrast 
to previous suggestions implying chemical bonding, we show that (02)4 is a van der Waals like 
cluster where exchange interactions preferentially stabilize the singlet state. However, as the clus- 
ter shrinks, there is an extra stabilization due to many-body interactions that yields a significant 
softening of the repulsive wall. We show that this short range behavior is a key issue for the 
understanding of the structure of e-oxygen. 
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The nature of the bonding in molecular oxygen clusters has been a subject of debate for 
nearly a century starting with the suggestion by Lewis of formation of dimers to explain 
the deviations from Curie's law observed in liquid oxygen[l|. A chemically bound dimer 
was expected considering the open shell character of O2 where two unpaired electrons 
occupy degenerate n* orbitals. However, a number of experimental {2 4 1 and theoretical jsj- 
Q] works clarified that in the gas phase (02)2 has the typical features of a van der Waals 
complex: a well depth of tens of meV and retention of the molecular properties within the 
complex. In fact it has been shown that a singlet species of D2h symmetry is stabilized 
due to exchange interactions but not in a sufficient extent to lead to chemical bonding. 
On the other hand, formation of a O4 molecule with four equivalent single bonds has been 
theoretically predictedjs] but roughly 4.3 eV above its (02)2 van der Waals counterpart. 
This O4 molecule is structurally closer to the naturally occurring sulfur rings Sg and in fact 
an analogous crown-shaped Og cluster has been predicted jol but again, as a very high energy 
isomer. 

Recently, the determination of the structure of the high pressure e phase of oxygen 10l.lll| 
has risen the interest in the study of molecular oxygen oligomers. In contrast to previ- 



ously proposed structures based on the dimer 
pendent x-ray diffraction experiments 1^, 



12| and herringbone chains [13], two inde- 



definitively concluded that e-02 consists of 



layers of well-defined (02)4 aggregates. They were found to form prisms with the O2 axes 
perpendicular [lo| or nearly perpendicular to a rhombic, nearly squared, base. There 
is a hierarchy of distances in this phase, the O2 bond length which nearly keeps the gas 
phase value (~ 1.21 A) at all pressures, and the intra- and inter-cluster distances (2.34 and 
2.66 A, at 11.4 GPa), which decrease monotonically with pressure up to the boundary with 
the metallic C-phase[ll|. Other key properties of this phase, suggesting increasing inter- 
molecular interactions, are a dark-red color, a strong infrared absorption and a magnetic 



collapse (12|, llJ-l^. Further evidence for a new intermolecular bonding has come from in- 
elastic x-ray scattering [17| where, at the lower pressure boundary of this phase (10 GPa), a 
discontinuous shift of about 1.1 eV in the electronic transitions from Is to Ivr* orbitals was 
found. 

A few works 
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18 



20| have attempted to rationalize the stability and the bonding of 



the (02)4 species in the framework of the density functional theory (DFT), but with unclear 
results. Thus, authors of Ref.jisl found that the D^h cuboid structure corresponds to a 
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local energy minimum that they recognized as unstable when higher levels of theory were 



applied. Furthermore, DFT calculations in Ref. 
(02)4 geometry is the most stable one compared with other chain structuresll3 
the need for additional studies. Then, it is apparent that despite recent progress [21 



19| failed to show that the experimental 

showing 
1|-|23| made 

in DFT methodologies for treating dispersion forces, the multiconfigurational character of 
molecular oxygen clusters still remains a serious challenge for such techniques. 

We report here high level supermolecular ab initio calculations of the (02)4 cluster. Our 
goal is a reliable characterization of the sin g l et g round state which is consistent with the 



magnetic collapse 



161 ] and spectroscopy 
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15| of the e phase. To this end, we use a mul- 



ticonfigurational ansatz which is unavoidable for spin multiplicities of (02)4 lower than the 
maximum one (nonet). We proceed in analogy to our previous work on the dimerjo], 7| and 
treat the highest spin complex by means of a restricted coupled cluster theory with singles, 
doubles and perturbative triple excitations [RCCSD(T)]. In addition, the singlet-nonet split- 
ting can be well described at the multiconfigurational complete active space second order 
perturbation (CASPT2) theory. Finally, the (02)4 singlet energy is obtained by adding, to 
the RCCSD(T) nonet potential, the singlet-nonet CASPT2 splitting. The aug-cc-pVQZ 24] 
basis set has been used at all levels of theory. For the CASPT2 calculations the active space 
is defined by distributing 8 electrons in 8 molecular orbitals correlating asymptotically with 
the O2 TT* shell. As customary, these orbitals have been previously optimized with the Com- 
plete Active Space Self Consistent Field (CASSCF) method. The counterpoise method {25 1 
was applied to correct interaction energies for the basis set superposition error. As for the 
(02)4 geometry, we use a cuboid structure with D4^h symmetry. The centers of mass of O2 
form a square whose side is changed in the range [1.5 — 25] A while the intramolecular dis- 
tance is kept fixed to 1.2065 A. Calculations have been performed with the Molpro2006.1 
package 261]. In addition and in order to study the role of many body interactions, we 
compare the supermolecular calculations just described with estimations based on the sum- 
mation of pure pair interactions. In the pairwise approach we have obtained expressions 
for the (02)4 energy which are compatible with a well-defined total spin of the complex, 
resulting in specific combinations of the (02)2 singlet, triplet and quintet potentials. For a 
faithful comparison, the pair potentials were obtained at the same level of theory than those 
of (02)4- It niust be noted that there are three singlet states{27| asymptotically correlating 
with four 02('^S~), and that here we are reporting the calculations of the ground singlet 
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FIG. 1: Upper panel: Interaction energies (in meV) for the nonet and ground singlet states of 
(02)4 as functions of the square side d (in A ). Supermolecular approach is represented in solid 
lines while the pairwise approach in dashed lines. Lower panel: many-body interaction energy for 
the singlet state, AV{d), obtained as the difference between supermolecular and pairwise energies. 
The large values of AV at small d's might be a key feature to explain the clustering of O2 molecules 
in the e phase. 



state. 

The interaction energies of (©2)4 in the singlet and nonet states as functions of the square 
side d are reported in the upper panel of Fig. 1 together with the pairwise estimations of the 
corresponding interactions. Equilibrium parameters of these potentials are given in Table 
I. As can be seen from the parameters of the singlet and nonet potential wells, (©2)4 is a 
van der Waals like complex in the gas phase, mainly stabilized by dispersion interactions. 
Exchange interaction, however, plays a role making the potential well of the singlet state 
deeper and shifted at shorter intermolecular distances than that of the nonet state. As in 

nn 

the dimerp, 17[, the exchange interaction favors the states of lowest spin multiplicity. 
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TABLE I: Binding energy D^, (in meV) and equilibrium distance (in A) of the (©2)4 potentials 
reported in Fig. 1 



nonet singlet 

De Re Df, Rf, 

supermolecular 56.6 3.37 122.6 2.99 

pairwise 57.1 3.37 108.4 3.08 



More insight is gained when comparing the supermolecular calculations with the pairwise 
estimations. For the nonet state, the pairwise approximation reproduces very well the 
supermolecular energies indicating that many-body effects are not particularly relevant. 
However, for the singlet state, an analogous agreement is only achieved for the larger d sizes 
of the cluster. Around the minimum of the singlet well, the supermolecular energies are 
already lower than the pairwise ones (see Table I) but it is for shorter distances where a 
remarkable softening of the repulsive wall is found. This is due to a many-body effect, as 
shown in the lower panel of Fig. 1 where it can be noticed that the many-body interaction 
energy for the singlet state increases dramatically as d decreases, beings about 1.2 eV at d= 
2 A. The origin of this effect must be in the exchange interactions [28;] since polarization 
contributions to many-body interaction energies, if important, should be shown also for the 
nonet state, that is not the case. 

In a molecular crystal the first response to pressure is the "squeezing out of van der 
Waals space" , in other words, the penetration to the repulsive region of the intermolecular 
potentials 29| . Therefore the peculiar behavior reported above can be relevant indeed for 
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For this reason, in the 



understanding the structure of the high pressure e phase 
following we will focus on the features of (62)4 at short cluster d sizes. 

Some clues for the softening of the repulsive wall of the ground singlet state are given 
in Fig. 2, where we report CASSCF calculations of the molecular orbitals arising from the 
interaction of the eight half-occupied vr* orbitals of O2. An analogous calculation for the 
singlet (02)2 is shown for the sake of comparison. The properties of these optimized or- 
bitals barely change from the asymptote up to the equilibrium distance {d ^ 3 A), but for 
shorter d^s the interaction does give rise to bonding and antibonding orbitals. Interestingly, 
the four bonding orbitals in (©2)4 are more stable than the bonding orbitals of (©2)2 and 
the associated occupation numbers increase faster as d decreases, as a result of many-body 
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FIG. 2: Upper panel: Energies of the HOMOs originated from the O2 vr* orbitals vs. intermolecular 
distance d, for the ground singlet states of (©2)4 and (02)2- Note that for the tetramer each line 
corresponds to a couple of almost degenerate HOMOs. Lower panel: As in the upper panel but 
for occupation numbers. 

exchange effects. However, a double occupancy only occurs for very short distances (d < 
1.8 A). Moreover, the (02)4 orbital stabilization is much smaller than the electron-electron 
Coulomb repulsion contribution to the total interaction. Thus, the result is not a mini- 
mum but a softening of the repulsive region of the potential, as obtained in Fig. 1. This 
multiconfigurational analysis differs from those of Refs. [l^, 3| where, based on a simpler 
monoconfigurational picture, it was suggested that all the bonding orbitals were doubly 
occupied leading to large binding energies with respect to the isolated O2 molecules. The 
present analysis also gives a qualitative insight into the observation of a 1 eV shift in 



the TT* ^ 1 s transitions at the boundary of the e phase [17|, since the splitting between the 
(02)4 antibonding energies and the isolated vr* orbitals is of the same order of magnitude 
1.5 eV) in the relevant range {d ~ 2.4 A). 

We have also checked that the energy of the ground singlet state is significantly lower than 
those of the other singlet states and different spin multiplicities as well, particularly in the 
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FIG. 3: Upper panel: Energy of a layer of e-02 (Eq. 1) as a function of the cluster size, d, for a 
pressure of 18.3 GPa (solid line). In the inset, the model unit cell is shown where the O2 axes are 
perpendicular to the a-b plane and the inter-cluster distances, rij, are displayed by dashed lines, 
being D the shortest one among them. A minimum of the energy is obtained for the intra-cluster 
distance dm- Lower panel: Pressure dependence of dm and Dm (lines) compared with data of 



Ref. 



Ill (circles) . See text for details. 



repulsive region of the interaction. Albeit in a very intuitive manner, the emerging picture 
can give some hints on the formation and structure of the e phase and its non-magnetic 
character. To illustrate this point and also in order to study whether present results are 
adequate for a description of the e phase, we consider a very simple model for the energy 
of a layer of (©2)4 clusters. A basic unit cell is shown in Figj3] where the tetramers form 
rhombuses of length d and angle a. It is assumed that, at a given pressure, the centers of 
mass of the clusters are fixed and determined by the lattice parameters a and b. Values of 
these parameters as functions of pressure were taken from Ref. In addition, the angle a 
is fixed to 81.4° as derived from data reported at 11.4 GPa ll|. We assume that all clusters 
increase/decrease their size at a time and study the subsequent modification of the cell 
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energy. This energy is given as a sum of intra- and inter-cluster contributions 

£;(d) = r"*™(d) + i5^y,f-(r,,), (1) 

where is the aheady reported supermolecular (02)4 singlet potential and V^*"**^^ is a 

pair potential between molecules i and j belonging to different clusters. The corresponding 
intermolecular distance rij is determined by d, a, h and a. A spin-averaged (02)2 potential 
was used for V^*"-*"^^ because the nonmagnetic character p^] of the e phase suggests that 
dependence on spin must be washed out. In Fig|3] (upper panel) it is shown that, within 
the cell, the optimum size of the cluster is considerably reduced with respect to the gas 
phase equilibrium distance. For 18.3 GPa, a minimum in the total energy is obtained at 
about dm= 2.17 A, in agreement with the observed intra-cluster distance of 2.185 A 11| as 
well as with the value obtained at 17.6 GPa in other independent experiment 10|]. In the 
lower panel of Fig. 3 we show that the pressure dependence of the optimized intra- and 
inter-cluster distances dm and Dm agrees fairly well with the observations [ill {Dm is the 
shortest inter-cluster distance, obtained as a function of dm-, a, b and a). 

We would like to stress that, despite the simplicity of the model, the key element is 
the behavior of the repulsive wall of the ground singlet state, adequately calculated at a 
multiconfigurational level of theory. Indeed, substitution of the the ground singlet energy 
ymtra^^^ with that Corresponding to a different spin multiplicity (or an excited singlet state) 



would lead to an optimum intracluster size far less compati 
dm ^ 2.35-2.40 A at 18.3 GPa). As discussed in Ref. 
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Die with the measurements (e.g., 
the fact that both intra- and 



inter-cluster distances compress at nearly the same rate is a clear indication of a rather 
weak interaction between the O2 molecules. Present finding of a van der Waals cluster with 
a very incipient chemical bond for short sizes is consistent with the observations. More 



refined models and extended ca. 
spectroscopy of e-oxygenjl2l. 14 



culations should be developed to account for the intriguing 



15| . Since at the cluster sizes relevant to the e phase the 



energy of the (©2)4 unit is not a minimum, including both intra and inter-cluster degrees of 
freedom in such models is unavoidable, especially for the lower frequency vibrational modes. 
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